Kalb-Ramond axion production in anisotropic string cosmologies 
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We compute the energy spectra for massless Kalb-Ramond axions in four-dimensional anisotropic 
string cosmological models. We show that, when integrated over directions, the four-dimensional 
anisotropic model leads to infra-red divergent spectra similar to the one found in the isotropic case. 
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I. INTRODUCTION 

The pre-big-bang (PBB) model of cosmology jl| in- 
spired by the duality properties of string theory, is faced, 
on the phenomenological side, with the question whether 
or not it can reproduce the amplitude and slope of the ob- 
served temperature anisotropy spectrum (|] and of large- 
scale density perturbations. 

Within the PBB model, the inflationary expansion due 
to the dilaton field in the low-energy effective action of 
string theory, leads to an amplification of metric fluctu- 
ations as well as of quantum fluctuations of all the fields 
present in PBB cosmology. Such fields, which are not 
part of the homogeneous background whose perturba- 
tions we study, are for example the gauge fields and the 
pseudo-scalar partner of the dilaton field in the string 
theory effective action. 

At first, it was thought that the PBB scenario 
could not lead to the observed scale-invariant Harrison- 
Zel'dovich spectrum of perturbations at large-scales. 
First-order scalar and tensor metric perturbations were 
found to lead to primordial spectra that grow with fre- 
quency jg], with a normalization imposed by the string 
cut-off at the shortest amplified scales. These blue spec- 
tra have too little power at scales relevant for the ob- 
served anisotropies in the cosmic microwave background 
(CMB). In contrast, the axion energy spectra were found 
to be diverging at large scales, red spectra, leading to 
very large CMB anisotropies, in conflict with observa- 
tions. 

These results already rule out four-dimensional 
isotropic PPB cosmology. However, if one allows for in- 
ternal contracting dimensions in addition to the three ex- 
panding ones, the situation is different. The axion field 
can lead to a flat Harrison-Zel'dovich spectrum of fluctua- 
tions for an appropriate relative evolution of the external 
and the compactified internal dimensions [HJH]. Thus, it 
is possible that the amplification of quantum fluctuations 
of fields which are present in the PBB scenario, can gen- 



erate via the seed mechanism || the observed anisotropy 
of the CMB radiation. 

Considering an isotropic PBB model with extra dimen- 
sions, the amplification of electromagnetic vacuum fluc- 
tuations and of Kalb-Ramond axion vacuum fluctuations 
lead to interesting observational consequences within the 
context of primordial magnetic fields [Q and large-scale 
temperature anisotropies In particular, massless ax- 
ions as well as very light axions can exhibit a flat or 
slightly tilted blue spectrum which may reasonably fit 
the observational data |,|. (Even though an acoustic 
peak at £ ~ 350 is excluded by experiments published af- 
ter Ref. H was completed, it is possible to shift this peak 
to I ~ 220 by closing the universe with a cosmological 
constant. More details about this model can be found in 
Ref. 0.) 

Recently it has been suggested that four-dimensional 
string cosmology models which expand anisotropically 
can also lead to blue or flat energy spectra for axionic 
perturbations [JTTfl . According to Ref. |i~lf , one can in- 
stead of assuming internal extra dimensions |gj , consider 
an anisotropic four-dimensional background. This has 
become especially interesting in view of new results which 
show that the pre-big-bang phase may generically be ho- 
mogeneous but anisotropic |12|. 

In Ref. [p!l| , the axion spectrum is only computed for 
the part of phase space where the longitudinal compo- 
nent of the wave vector is sufficiently large. In this work 
we correct the result of Ref. |ll| and complete the com- 
putation to contain all directions in phase space. We 
then integrate the obtained spectrum over directions and 
compare it with the result for the isotropic PBB. We find 
that the anisotropic spectrum, when averaged over direc- 
tions agrees roughly with the isotropic one. Therefore, 
anisotropic expansion during the pre-big-bang phase can- 
not solve the axion problem of four-dimensional string 
cosmology. 
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II. AXION PRODUCTION IN THE 
P RE-BIG-BANG COSMOLOGICAL MODEL 

Let us consider a four-dimensional spatially flat 
anisotropic PBB cosmological model, with metric 



( 5 ^) = diag[l,-a 2 (t) ) -6 2 (t),-6 2 (t)] 



(1) 



the internal compactified radii (if present) are assumed 
to be frozen. For simplicity, we assume two directions to 
expand with the same scale factor b. Varying the low- 
energy string theory effective action (in the string frame) 
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(H^ va denotes the antisymmetric tensor field), with re- 
spect to the metric and the dilaton field cj>, we obtain 
the dilaton driven vacuum solutions of the tree-level evo- 
lution equations. As derived in Ref. p|, these solutions 
read 



H 
L m 



and 



<P(v) 



a + 2(3 - 1 



log 



m 



1 — a 

with a and (3 satisfying the Kasner condition 
a 2 + 2(3 2 = 1 . 



(3) 



(4) 
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Here rj denotes conformal time with respect to the scale 
factor a. It is negative during the pre-big-bang era and 
77 = —rii stands for the transition time from the dilaton 
driven pre-big-bang era to the radiation dominated post- 
big-bang era. To obtain the axion evolution equation, 
we vary the effective action, Eq. (0), with respect to the 
Kalb-Ramond axion field a, given by 



-d a a . 



(6) 



The evolution equation of the canonical field ip — e^/ 2 ba, 
in Fourier space, reads (ll]] 



k 2 — 



i>k = o 



^l 2 b , (7) 



where kh denotes the modulus of the comoving longitu- 
dinal momentum and kx — W 1 ky + k^ is the modulus of 

the transverse momentum. Equation (g) describes the 
generation of axionic modes, where the anisotropy of the 
spacetime has been translated into an asymmetry be- 
tween the longitudinal and transverse momenta. 



The choice a = (3 — — l/v3 corresponds to the 
isotropic case, for which a — b and the evolution of ax- 
ionic fluctuations is given by 



Vk 



with ij) k = Va k , V 



ipk = , 
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. V" (/J 2 -1/4) . . 2 , 1 N2 
so that — = ^ , with fi = (p — -f 



(8) 
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The solution of Eq. (g) , normalized to an initial vacuum 
fluctuation spectrum, can be written as 



1 p k = r) 1/2 Hj?H\W) > A* 



, V < ~Vi , (11) 



with p = Hjj, denotes the Hankel function of second 
kind (we adopt the conventions of Ref. |l3| ) . 

Assuming that the dilaton driven era is followed by a 
radiation dominated era, the density parameter of pro- 
duced Kalb-Ramond axions per logarithmic frequency in- 
terval is [0 



0<7 (w, 77) 
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1 dp a 
p c d log uj 
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where p(u>) denotes their spectral energy density and 
p c = 3M 2 H 2 1 (8ir) stands for the critical energy den- 
sity. Note that wi — k\/a\ — l/(a 1 |?7 1 |) represents the 
maximal amplified frequency, g\ = H±/M p is the transi- 
tion scale in units of the Planck mass, Hi ~ uj\ denotes 
the Hubble scale at which the universe becomes radia- 
tion dominated. Hence ^7(77) = (Hi / H) 2 (a\ / a) 4 is the 
radiation density parameter at a given time 77. 

Clearly a flat spectrum corresponds to p = 3/2 and the 
value p = v3 obtained in a four-dimensional isotropic 
pre-big-bang model implies a red spectrum, leading to 
an unacceptable divergence at low frequencies. 

Let us now go back to the case of a four-dimensional 
anisotropic background. We first study the evolu- 
tion of axionic fluctuations and we then calculate the 
spectral energy density of the axionic inhomogencitics 
(dp a I dlogu)), as they re-enter the horizon during the 
isotropic radiation dominated era, after being ampli- 
fied during the anisotropic dilaton driven era. Inserting 
Eqs. (g) and (g) into Eq. (g), we obtain jD) 



< + [kl + k\ 



771 



- 1/4 
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where 



2(a-,3) 
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If 7 < 0, the kx-teim as well as the 77 _2 -term go to zero 
for j] — > — oo; and initially the parentheses in Eq. is 
dominated by k\ (except if = 0). If kr is not very 
large, namely if 



k T < Mfci/M~ 7/2 



(16) 



the scale fct becomes super- horizon, i.e. the parentheses 
in Eq. ( |l3| ) is dominated by the l/?7 2 -term, before the 
fcr-term takes over. In this case, we may entirely neglect 
the fcr-term in Eq. (|l3|), which then reduces to Eq. (||) 
with k replaced by k^. Therefore, the spectrum for these 
modes is flat for /i = 3/2 which corresponds to 

a = -7/9 , /3 = -4/9 and 7 = -3/8 . (17) 

We also require the solution to expand, i.e. a, (3 < 0. 
We first concentrate mainly on these values of the Kasner 
exponents since they lead to a scale invariant spectrum 
of fluctuations for directions with a sufficiently large k^- 
component, but we express our results in terms of a and /3 
so that they can then also be applied also to other values 
of the Kasner indices. In the part of phase-space defined 
by the inequality given in Eq. (^), the energy density 
of the produced axions has already been determined in 
Ref. fll]| . Here we correct the result of Ref. 11 and 
generalize it to the entire phase space. 

To solve Eq. (|ll]), we distinguish among the following 
two cases: 

(I) The modulus of the longitudinal momentum, fc^, al- 
ways dominates until rj 2 < 1/fc 2 at which point the l/rj 2 
term comes to dominate. This is equivalent to the con- 
dition given in Eq. (16). 

(II) At some conformal time r\ = t\t < — 771 , the modu- 
lus of the transverse momentum, kx, comes to dominate 
over fc^, but the mode is still well within the horizon, i.e. 
\/k 2 L + k\{—r)T Mi) 1 > Vt 2 - Equation ( |l3| ) implies 



Vt 



'k L 
11 { .h 
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Case (I) : Let us first discuss this case which is also 
the one studied in Ref. jyj. Here, the inequality given 
in Eq. ( |l6| ) holds. For low frequency modes, u> <C u>\ 
this is the case outside a very thin slice around the plane 
fci = 0if7<0. In this situation we may entirely neglect 
the second term inside the parentheses of Eq. ( J13| ) which 
yields a Besscl differential equation. Its solution during 
the pre-big-bang era , is simply 



for 77 < —Tlx 



(19) 



After the transition to the radiation dominated era 
(RD), we assume the dilaton to be frozen and the ex- 
pansion to have become isotropic. This implies V" = 0, 
a/b = 1 and Eq. (Q) reduces to a simple harmonic equa- 
tion with general solution 



1 



c+e - ifc (»i+''7i) _|_ c _ e ik(v+m) 

for 77 > —T]i . 



(20) 
(21) 



By matching the in-coming solution i/^ BB to the out- 
going one V'fe ? an d by also matching their first deriva- 
tives, at the transition time 77 = —771, we obtain the fre- 
quency mixing coefficient c_(fc): 




{kT)i)(k L rn)^ 



(22) 



The coefficient c_ determines the occupation numbers 
of produced axions. The spectral energy density of the 
produced axions reads 



Pl(oj,s) 



dp a 



-Mw)| s 



(23) 



FromEqs. (22) 



PL (w,s) 



d log U) 7T 2 

) we obtain with oj = k/a for /j = 3/2 
1 



2tt 3 



(24) 



where s — k^/k. 

Thus, if the longitudinal momentum k^ dominates, the 
spectrum of produced Kalb-Ramond axions is flat, i.e. 
independent of u>, but anisotropic. This result generically 
agrees with the finding of Ref. [jlljl (up to a factor 1/s 2 , 
which we think is missing in Refi ll l[| ) . 

Case (II) : We now assume that the kT-teim comes 
to dominate before the perturbation becomes superhori- 
zon. As long as the perturbation is sub-horizon, we may 
approximate Eq. (O) by 



m 



1 T 



i> k = 



(25) 



(18) An approximate solution to this equation is 



exp \ri\Jk\ + q 2 {—rj/rii)~ik 



^[kl + i-v/TIi)^ 4 



(26) 



with q = 1/(1 + 7/2) = (l-a)/(l-/3). 



In the regime considered, r]\Jk 2 L + q 2 {— ?7/?7i) 7 fcy 3> 1, 
this solution becomes exact, if either fc^ or kr vanishes 
and it is a good approximation if one of the two terms 
dominates. If the fc^-term and the fcr-term are of the 
same order, the relative error is about |7/2| = 3/16. It 
is also clear that this represents the correctly normalized 
incoming vacuum solution. 

At conformal time 77 = tjt, the transverse momentum 
kx comes to dominate over the kjy term in Eq. (|T3|). At 
even later times, the ?y _2 -term will eventually dominate. 
After r/T Eq. ( |13[) can be approximated by 









. Vi. 



7 ll 2 



1/4 



(27) 



3 



with general solution 13 
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where q is as above, and 4^ are Hankel functions 
of the 1st and 2nd kind of order fiq. For large kr\r]\ the 
second term just corresponds to the solution (|2^) in the 
limit where fcz, can be neglected. Therefore, by matching 
the solutions we find 



c« = 

J 2 ) _ 



up to an irrelevant phase. 

Next, we have to match the solution for the field a 
of the pre-big-bang era to the solution for a during the 
radiation era at the transition time r\ = —Tji. 

As we go from the pre- to the post-big-bang era, we 
assume the universe to become isotropic and the dilaton 
field 4> to become frozen. Thus, here the matching of the 
in-coming to the out-going solution for a, is not equiva- 
lent to matching ip. The relation between the the axion 
field a and the canonical field ip at conformal time r\ is 



(30) 
(31) 



The canonical field in Fourier space during RD is given 
in Eq. (|2^) . Matching the solutions and their first deriva- 
tives for a, as we pass from PBB to RD at time rj = —Tji, 
we obtain for |fori7i| "C 1, the Bogoliubov coefficient c_ 
given by 



-2fj,q 



XS 



-2p,q 



(32) 



Inserting the values /i = 3/2, a = — 7/9 (3 = —4/9 which 
lead to a flat spectrum in case I one finds a somewhat 
blue spectrum in case II, 



Pt(u, s) (xuj 9/13 . 



(34) 



(28) Of course this case also gives a finite answer on the 
plane kh = for which the result obtained under case I 
diverges. 



(29) p(u,s) 



III. RESULTS AND CONCLUSION 

In total we can summarize the calculated spectrum by 



-7/2 



2tt 3 



3-2/xg 

(l- s ^-m ^ else, 



(35) 



where s = k^/k , p, = 1 — 2/3/(1 — a) and q — (1 — a)/(l — 
j3). For our prefered values, a — —7/9 and j3 = —4/9 
which imply [i = 3/2 and q = 16/13, the above result 
reduces to 



s- 3 if k T <k L (k 1 /k L )-"</ 2 

■>a / \ -9/13 

(Is 2 )-^ (%) else. 

(36) 



In the regime of phase space where the longitudinal mode 
of the momentum is very small i.e. when the condition 
given in Eq. ( |l6| ) is violated, the spectrum of the pro- 
duced Kalb-Ramond axions is not flat. For a given value 
of u>, this is the case if s is smaller than the critical value 
s c which is well approximated by 



s c (w) ~ — 



3/13 



if uo ~ O.lwi 



(37) 



a very small value for cosmologically interesting frequen- 
cies. 

In Figure 1 the the energy density p(u>) is shown as a 
function of s for different values of uj. 



With Eq. (|23j) we then obtain that the energy density 
of the produced Kalb-Ramond axions, in the case where 
the transverse momentum kx comes to dominate, i. e the 
inequality given in Eq. ( |l6| ) is violated: 



Pt(w,s) 
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(33) 



4 



-10 



-6 



-4 



-2 



l°g 10 s 



-30 



-20 



10 



FIG. 1. The density parameter of the produced 
Kalb-Ramond axions is shown in units of gl0. 1 as a function 
of s for different values of w. The value uj\ ~ (0.01 — 0.1)Mp; 
is the string scale. For this curve the results obtained in 
Eq. ( |36[ ) have been interpolated logarithmically. The little 
bars on the log(s)-axis indicate the values of s c at which the 
condition (^) becomes an equality. For values of s smaller 
than Sc(ui), the spectral density is approximated by pt(ui). 

For s fixed, if the modulus of the longitudinal momen- 
tum dominates in Eq. (|ll|), more precisely if it satisfies 
the condition fc^ > fc*r(&TAi)~ 7 / 2 ^ 1+7 / 2 \ the spectrum 
of the produced Kalb-Ramond axions is flat. 

To estimate the total energy density per logarithmic 
frequency interval we integrate the axion density Sl a (uj, s) 
over s. For this we use 



d 3 k = 2nk T dk L A dk T = Airk 2 ds A dk , 
where we have used dk^ — kds + sdk and 

kds + \J 1 — s 2 dk . 



(38) 



dkx = 
Hence, we have 



Q, CT (u),r)) = / Q a (ui,s,ri)ds 



1 

Pc 



Pt(w, s)ds 



Pl (w, s)ds 



/Wl\9/13 0.5 

V uj J si 



(39) 
(40) 



This spectrum is shown in Fig. 2. 



log(u/"i) 



FIG. 2. The energy density D. cr (ijj)/g 1 D. 1 , integrated over 
directions s is shown as a function of uj (solid line). Comparing 
it with the isotropic result (dashed line) we conclude that, the 
two spectral indices are the same and, within our accuracy, 
also the amplitudes are comparable. 

Using s c ~ (w/wi) 3 / 13 , which is a good approximation 
as long as uj < 0.1 wi it can also be seen directly from 
Eq. (M) that the isotropic energy spectrum is nearly 
reproduced. The isotropic spectral index, 3 — 2\/3 ~ 
—0.464 is actually replaced by —6/13 ~ —0.463. Insert- 
ing reasonable values for the string scale, 0.01 < g\ < 1, 
we see that also in the anisotropic case axions are over- 
produced in unacceptable amounts. Even if the spectrum 
of the axions from wave vectors directed sufficiently far 
from the plane k^ — 0, is scale-invariant, the enhance- 
ment of the spectrum in the vicinity of the plane k^ = 
leads to a total contribution which agrees with the one 
obtained in the isotropic case. Therefore, the model is 
excluded (see Ref. 0). 

So far we have mainly considered the case a — —7/9 
and f3 = —4/9, but our results apply quite generically, 
as long as 7 < and thus the fc^-term dominates at 
sufficiently early times. But also if 7 > 0, Eq. ( p^ ) is 
an approximate solution on sub-horizon scales. In this 
situation, however the fcr-term dominates at sufficiently 
early times and continues to do so until the perturbation 
becomes super-horizon if the inequality given in Eq. ( |l6| ) 
is violated. For 7 > this is the case outside a narrow 
cylinder around the kx = axis. Therefore, the generic 
formula given in Eq. ([39]) always applies, but s c <C 1, if 
7 < and s c ~ 1, if 7 > 0. 

For general values of a and (3 we obtain 



f2 (r (w,r?)= / tt a (uj,s,r])ds 



(1 - s 2 y^ds 
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3-2/x 



s c (w) 



The transition value of s is given by 



7/2 



(41) 



(42) 



/ \7/2 

If 7 < (i.e. a < p), the factor ( ^ 1 is very large 

in most of phase space and hence s c < 1. On the other 
hand, if 7 > (i.e. a > (3), the above factor is very 
small for the relevant frequencies, uj <C lo\ and s c ~ 1. A 
reasonable approximation is 



q-l 



1 



2/9-2 



if 7 < 
if 7 > 



(43) 
(44) 



where we have used the relation q = 1/(1 + 7/2). Insert- 
ing these results in Eq. (J4l|), the integrals can be approx- 
imated by 



where 



(45) 



n =2 + q- 2 f iq = 1 + a + 2 P if a < p _ (46) 
1 — P 

n = 1 + 2/g - 2 M = 1 ± - + 2/3 iia>0. (47) 
1 — a 

Clearly, since a 2 + 2/3 2 = 1 and a, (3 < it is a + 2/3 < -1. 
This shows that the spectrum is never blue and becomes 
scale invariant only in the degenerate case with two static 
dimensions, (3 = 0. This is also shown in Fig. 3, where 
the above approximation for the spectral index plotted as 
a function of a: the spectrum always remains red with 
a spectral index relatively close to the isotropic value, 

^iso — 

3 - 2\/3 0.46, except in the extremal case, 

when two dimensions are frozen and a = — 1. 

If one relaxes the condition that both a and b be ex- 
panding and just asks for volume expansion, a + 2(3 < 0, 
there is another pair of values for the Kasner indices lead- 
ing to a flat spectrum, namely a = —1/3 and (3 = —2/3. 
However, if we want expansion in all three dimensions 
the spectrum is always red. 




FIG. 3. The spectral index n is shown as a function of 
the exponent a of the expansion law. For a = — 1/V3, the 
isotropic case, our approximation is not very good since there 
a — j3. This is reflected in the unphysical kink at this value 
of a. Clearly, the resulting spectrum is always red (n < 0), 
with —0.4 > n > —0.5 except close to the degenerate case 
a -> -1. 

To summarize, we find that anisotropic expansion has 
very little influence on the overall axion production and 
cannot cure the axion problem of four-dimensional pre- 
big-bang models. Only by allowing for extra dimensions 
one can escape this conclusion and obtain a scale invari- 
ant spectrum of axions as described in Refs. A 
'realistic' string cosmology with a Kalb-Ramond axion 
can therefore be realized only in models with extra di- 
mensions. 
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